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SUMMARY 


This paper presents a study of the nature of the 2k hour 
synchronous satellite perturbations due to the earth’s triaxiality. 
Equations representing the drift due to these perturbations are 
developed. It is found that the radial drift is linear and the 
longitudinal drift is parabolic. Both have superimposed on these 
primary curves, oscillations with one day periods. This analysis 
also verifies the existence of only two dynamically stable points, 
where no drift due to the earth’s triaxiality will occur, 
are located at opposite ends of the equatorial minor axis. 
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INTRODUCTION 


It has been^shown (Refs. 1 and 2) that there are only two dynamically stable 
points for 2k hopr synchronous satellites. These sure in the vicinity of 57° east 
longitude and 123° west longitude in the geographic system. Satellites posi- 
tioned at locations other than these will drift toward the nearest point of 
stability. In order to obtain maximum global coverage using the synchronous 
communications satellite concept, it is necessary to maintain at least three 
satellite systems, which sure equally spaced about the equator, in operation. 

This means that drift correction is necessary. General perturbation studies 
have shown the disturbances of satellite orbits to be caused primarily by the 
potential field of a triaxial earth. The indicated influence of the lunar and 
solar gravitational fields is small in comparison. It is the purpose of this 
paper to determine the nature of the 2k hour synchronous satellite perturbations 
due to the earth's trisociality. 


DEVELOPMENT 


REFERENCE SYSTEM AND DEFINITIONS 


An inertial reference system with axes X, Y, and Z as shown in Figure 1 
is used in this paper. The origin and Z axis of the coordinate system coincide 
with the earth ' s center and polar axis respectively. The X-Y plane coincides 
with the equatorial plane and completes a right hand coordinate system. It is 
found convenient to specify the satellite position by the spherical coordinates 
0, P and h. The earth's equatorial minor axis lies at an angle, 

©g- =* 6c + ©gt (l) 

o 

from the X axis, where 

© g = the initial angle between the X axis and the equatorial minor 
° axis 
• ' *- 

© £ = the earth's rate of rotation about its polar axis 

and, t = time 

From Figure 2 the relationship of inertial to geographic longitude is : 

© ='«£+ A M (2) 

where X = geographic longitude 

and, yQ = the angle between Greenwich (zero longitude) and the equatorial 

minor axis and is approximately 123 ° west longitude 

Further the angle between the equatorial minor axis and the projection of the 
satellite radius vector, r, into the equatorial plane is given by 

a* = a i /3 (3) 
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EQUATIONS OF MOTION 


The components of acceleration in terms of r, 9, and sere: 

Q. h » r - r©*cos 2 ^ - r</) 2 (*0 

«s * r cos ft M r ** cos M (5) 

* *r + rG^cos 0 S//7 <f> ( 6 ) 

If one assumes the influence of the lunar and solar gravitation fields to 
be small, the force components are determined from the earth's gravitational 
potential. The most recent formulation of this potential is: 


u 


_ J ±_ __ Jzo /< ftp 




3SI 




where 


zdzzM^Ro C oS*</> COS 

1^3 * 



^ = Universal gravitational constant 

/? 0 = Mean equatorial radius of the earth 
and v^20 an< ^ ^ 22. 8X6 ^i™®* 18 ! 011 !® 88 constants whose values are (Reference 2): 
J 20 - 1.082 X 10“ 3 
J 2 2 = “5*35 X 10“ 6 


It is seen that Eq. ( 7 ) includes one each of the zonal and sectarial harmonics. 
All higher order harmonics are assumed negligible. 

From Eq. (7) the force components are: 


F * * 


m ZU _ 
m JF - 


F a = 


+ ^2 £ \ ffi— ( ^ sirf^ - l) 

aJ^HRc, COS*ft COS 2^ 

r fe “ m (~ cos V s//9 ( 9 ) 


(e) 


2 



and. 



_ / so . 

m T Sf - 



sMsgfL sm<f> cos 


(10) 


where 

tn = Mass of the satellite 

By equating the corresponding acceleration and force expressions from 
Eqs. (4), (5), (6 ), (8 ), (9) and (10) the following equations of motion are 
obtained. 


2 \ 

— cosj) sih<t> cos 2$ 


r - re* cos* <f> -r<f>* = -^1 

- ij - aj&gdk. cos 2 <f> cos 2 $ (11) 

r ~cos<f> cos ^) = “ cofysin'2# (12) 

+ r£cos<f)sin<f> - - s/n <j> cos<j> 


— & 

h+ 


cos $ s/n <f> cos 2 y (13) 


General perturbation studies have indicated the latitude perturbations to be 
periodic and having a magnitude of less than one degree. For purposes of this 
report, latitude is assumed equal to zero at all times. Therefore, Eq. (13) 
is eliminated and the analysis becomes two dimensional. Eqs. (ll) and (12) 
reduce to: 


•• 'a _ _ _ 

r - re - ^ 


3^/20^ & 


2 M 


9fa42 /£ 
h+ 


COS 2# (l4) 


+- 2A0 


<S t/22 ^ /?o 

pr— 


S/h2t 
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(15) 



Since the primary concern is with <fewations of a satellite from a synchronous 
orbit, it becomes advantageous to linearize Eqs. (l4) and (15) in the vicinity 
of this orbit. Perturbations about the synchronous orbit will be radial, > 
and longitudinal, Wow define: 


= e e +■ V 0 + 

( 16 ) 

» ft + A*- 

(17) 

= r s +• At- 

( 18 ) 


where = the desired longitude difference between the equatorial minor 
axis and the satellite position 

and, - synchronous radius 

Substitution of Eqs. (l6), (17) and (l8) into Eqs. (l4) and (15) 

I 2 ... 


and division by 

r s ®e 

4/ 

At\- 






_ 3 Jio M Ro ( 1 - 4 £!- \ 

2 6/ V l r s) 




and. 


- 4- -^){cos2l - 2 At sm2i,J ( 19 ) 

+ 2 

= - - 5 [sinzt + 2 AX cos 2 ^ 1 . 20 ) 


where small angle approximations have been made and terms of second order and 
greater have been neglected. A dimensionless time is now introduced as: 


r = e £ t 


_ k - 


( 21 ) 


In terms of this dimensionless time, Eqs. (19) and ( 20 ) reduce to: 


jL 2 / 

■f) ' C (^ r )~ Z n A * ' * * 

(22) 

2 " F ft) + + EA * = 8 

(23) 

where 



A * 

1 - -4^3 - _ SJz.z'^fo cos 2 $ 

eFij? zefr? &i V 

( 24 ) 

B = 


( 25 ) 

C = 

1+ s?r 5 3 + e#r 5 ^ + #r/ co * 2a 

(26) 

D = 

iS x¥^ sm 2 & 

( 27 ) 

E * 

12 cos 2S ° 

(28) 

and P = 

— 2 & 

(29) 


By using operater notation, the variables of Eqs. (22) and (23) can be separated. 
The result is: 



(S 4 + fa s 2 +/^ 2 S +/S ) = AE + BD 

j 1 s 

( 30 ) 

and. 

(S 4 + ^ S 2 +/^ 2 S + 03 ) = AF - BC 

( 31 ) 

where 

S » d 



dr 

( 32 ) 


fa m 4 - C + E 



^ « 2 (D - F) 

( 33 ) 

and. 

@2 m - CE - DF 

( 34 ) 
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Examination of Eqs. (2^) thrc/ (29) shows the following to be very good approx- 
imations . 


A =0 

(35) 

b = -gJ 2Z sm 

(36) 

C =* 8.0 

(37) 

D = ieJ 22 . \^pr | ^ //7 2 #o 

(38) 

E = 12 Jzz |Aj*C0S 2^0 

(39) 

and, F ^ 30 J 22 ^Aj Sm 2ft 

(^0) 

Further , it is seen that: 


fit * 

(41) 

|S 2 s 24- /r« W7 2& 

(42) 

(S3 s 36/^ COS 2*. 

(43) 

/(£ + 8£> S 0 

(44) 

flF - SC = -18K^ SwZVo 

(45) 

where K 22 =-J £2 


Eqs. (30) and (31) then become: 

[s* * S z * (24 ^ 22 s/n2y o j $ 

+ (36 ^22 cos 2tf 0 j 

m 

Ar _ 

-p = O (46) 
'5 

+ s z + [24 K^smz^s 

4 

“4* A 22, 2-$>J 

Afr 

J . (47) 


- - is K zz sm z# 0 
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It is seen that the characteristic equation for both A* and -p* is : 


A = $ * + S Z + S -h j8 3 ( 48 ) 

Extensive numerical analysis of Eq. (48) for the values of and ^5 

involved, has shown 

A = [S z -t- p z 5 + /S 5 ) (s 2 - &S +i) (1*9) 

to be the solution as a product of two quadratics. The individual roots are: 

S t = lZK zz sm2» c + j (50) 

$ 2 =■ 12 K zz smzlf, -j (5i) 

S 3 = -J2^ i2 sm2 y o -f-J^ai^cosWo (52) 

S* = -12K 2Z Sln2& -J-3aK zz COS2% (53) 

For 0° 6 ^ < 45°, the characteristic solution is: 

A = e\ lzH2 - 2Sm2 *°\ z (c 1 stn? C z cost} 

+ C 3 e (~ 1Z *22 SM2X 0 +F 3G K^COSZf^ 

+ e (-*2 smzto cos 2 &)? (54) 

For 45° is 30^ the characteristic solution is: 

A = e^ 12/<zzSmZ ^ T ( i Ci s/nz + c z cosZ ) 
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+- e ("^ s//?2X 0 )7^c^ smi-z&K?.* cos 2 # 0 'Z 


+ C+ COS /- 36 /r 22 COS 2 &' 'T 


( 55 ) 


The exponential terms of Eqs. (5^-) and (55) have very large time constants 
and can, therefore, be replaced by unity. Also, the transcendental functions 
corresponding to the CU and C> coefficients of Eq. (55) involve extremely 
small angles and can be replaced by small angle approximations* The solution 
for all values of then becomes: 

A - C,$tnZ + C 2 cosZ + C z Z + C+ (56) 

In order to evaluate the unknown coefficients , the initial conditions must be 
determined. It is assumed that a perfect injection into a synchronous orbit 
has been achieved. This implies: 

Al r o = Al ~ = AV. = A& = O (57) 

Substitution of Eq. (57) into Eqs. (22) and (23) yields: 


A\Q * Atf 0 = O 


(58) 


'2~ z A){‘ 0 < 0 K zz sm2)f 0 (59) 

and, ^^,3^ ~ 1 2 K zz SIH (60) 

Evaluation of Eq. (56) and its derivatives for the 
for the characteristic solutions: 

(Atf) c = 6 K 2i S/n2*„[i - COST) 

and, = 12 sinZ&| Z - SlflZ) 

There is no particular (steady state) solution for 
For A* , the particular solution is: 

(Att) p = - (9K zz sm2&) Z z 


initial conditions gives 

(61) 

( 62 ) 

r* # 
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finally, the total solutions for Atf and AiT' , respectively, are: 

Ay = ^ • 6 k zz sm 2 y„ ( i - cos 4^ -l-segt*) (64) 
Ar = •V i2A iz- sv ' ,2 *«( 4* - (65) 

where real time has been substituted. As shown, the unit of A* is degrees 
and the unit of AT corresponds to that of . 


RESULTS AND DISCUSSION 

Ibr a truly synchronous condition to exist, the forcing functions of the 
coupled perturbation equations, Eqs. (30) and (3l)> must be zero. It follows 
immediately that A and B given by Eqs. (24) and (25), which are repeated below 
as Eqs. (66) and (67), must be zero. 


A » 1 - 


M. 3J20 M Ro 9^22 .^cRo s> r\e Vt 

zr?- rg c 1 — COa 2: <J 0 


2& 0 


rf 


elrrl 


( 66 ) 


B 


= - 6 sm2t a 


ri 


(67) 


Obviously, for Eq. (67) to be zero, sin 2$ 0 must also be zero. This gives 
four values of y o (0°, 90° > l80° and 270°) at which a truly synchronous condi- 
tion exists. Blitzer (Ref. 3) has shown that the locations - 0° and l80° 
represent dynamically stable points while the other two points are statically 
stable. The value of r which makes Eq. (66) zero at the dynamically stable 
points is found to be: 


= 22752 . 2 .S 2 - nautical miles (68) 

Shown as Figures 3 and 4 are plots of Atf and A\T , Eqs. (64) and (65), 
as functions of time for various values of ^ . It is seen that the satellite 
motion is symmetrical about & = 45°. For the cases shown, the drift was also 
determined by numerical integration of Eqs. (ll), (12) and (13 ) . Agreement of 
the two methods was very good within the limits of the linearizations used in 
the perturbation equations. These limits are: 

Ar << r s 

a* < 5° 
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(69) 

(70) 



The limit on A* is established by the limits of the small angle approxi- 
mations which were used. Figure 5 is a plot of the differences between the two 
methods of solution for = 45° where the drift magnitudes are maximum. It 
should be realized that for longer drift periods Eqs. (64) and (65) represent 
only short time period drift . Pbr example at values of r 0 in the vicinity of 
90° there is very little sh$^/|dr?f£° But since = 90° is only statically 
stable, it follows that long time period drift will be present. 

It is noted that if the oscillations are ignored, Eqs. (64) and (65) are 
identical to equations developed by Frick and Garber (Reference 2). The period 
of the oscillation in both cases is: 

PERIOD = % = * (71) 

When && is expressed in degrees, the maximum amplitude oi its oscillating term 
jLs * ^ 

AMPLITUDE = 6Jz2.[% “ .4208 X (72) 

The amplitude of the oscillation in degrees is very small. However, in terms 
of linear position this amplitude becomes 

AMPLITUDE = r s -6j zz [-&) = lOl. S ft. (73) 

The maximum amplitude of the AIT oscillating term is: 

AMPLITUDE = r s • 12 = 203 ft (74) 

It therefore seems that the oscillation terms are of sufficient magnitude to 
warrant consideration in drift correction studies. 

Figure 6 shows the nature of the drift in each of the four geographic 
quadrants. The symmetry is readily apparent. It is pointed out that the 
magnitudes of drift are equal in all quadrants. Only the directions are different. 

CONCLUSIONS 


The purpose of this paper was achieved in that equations were developed 
which accurately determine drift for short time periods. Consequently the 
following conclusions can be stated. 

1. The only locations at which truly synchronous conditions exist are 

in the vicinities of )f 0 = 0° and l80° at an orbit radius of 22752.292 

nautical miles. 

2. A realization of the presence of oscillatory motion in both the 
longitudinal and radial directions was achieved. This oscillation is 

of sufficient magnitude to warrant consideration in drift correction studies. 
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3 . The symmetrical nature of the drift in the four geographic quadrants 
is such as to greatly reduce the analysis of different longitudes as 
potential satellite locations. 
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Figure 3. Longitudinal drift 
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Figure 4 . Rrdirl drift 



»/> 

Ui 



1 

Figure 5 . Comparison of the drift results 
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